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A nonlinear susceptibilities (the third derivative of a magnetization ms by a 

> 

'nI" \ magnetic field /i ) of the 3=1/2 ferromagnetic Heisenberg chain and the classical 



Heisenberg chain are calculated at low temperatures T. In both chains the nonlinear 



■^ \ susceptibilities diverge as T ^ and a linear susceptibilities diverge as T ^. The arbi- 



B 

P^ ' function i<'(x)=(2x/3)-(44x'^/135) + O(x^) is common to all values of spin S. 

o' 



trary spin S Heisenberg ferromagnet [ 7i = J2i=i{~JSiSi^i — {h/S)Sf} (J > 0), ] 
has a scaling relation between nis, h and T : ms{T,h) = F{S'^Jh/T'^). The scaling 
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I. INTRODUCTION 

In the 1970's integral equations have been proposed to calculate thermodynamic quan- 
tities of the S=l/2 Heisenberg chain |l|J^. In this calculation the Bethe Ansatz method 
was used to diagonalize the Heisenberg Hamiltonian. One of the authors has calculated the 
Gibbs free energy numerically by the integral equations [^. 

On the other hand, the other approach has been developed. In this approach a quantum 
one-dimensional system is transformed to a classical two-dimensional system by the Trotter 
formula. Betsuyaku has defined a transfer matrix and he has calculated an energy and a 
specific heat numerically of the S=l/2 XY model and the S=l/2 Heisenberg model @]. By 
diagonalizing the transfer matrix with the Bethe Ansatz method, Koma f^ has calculated 
the Gibbs free energy and a magnetic susceptibility of the 5=1/2 Heisenberg chain. These 
quantities coincide numerically with results of ref.3. Moreover a correlation length of the 
S=l/2 Heisenberg chain has been calculated numerically [0. Furthermore Suzuki, Akutsu 
and Wadati pointed out that the transfer matrix of the 5"= 1/2 XYZ model was a special 
case of the inhomogeneous eight- vertex model which was solved by Baxter ||^ . One of the 
authors |p has obtained the Gibbs free energy and the correlation length of the S'=l/2 XXZ 
model in a magnetic field by solving self-consistent equations for infinite numbers, which 
were derived by the Bethe Ansatz method. 

In this paper we try to calculate the free energy of the 5*= 1/2 ferromagnetic Heisenberg 
chain in a magnetic field by solving self-consistent equations for infinite numbers. We also 
calculate the linear and the third order susceptibilities of the classical Heisenberg chain 
analytically. By comparing these results, we propose a scaling function of the arbitrary spin 
S Heisenberg chain defined by 

N f^ N 

n = -JJ2s^s,+,--Y.s^, (1.1) 

1=1 "^ 1=1 

where Si = {Sf,Sf,Si) are spin operators and J is the positive coupling constant. A 
partition function Zj^^Sy the Gibbs free energy per site fs and a magnetization ms of this 
model are defined as follows: 



Z^,5 = Tr|exp(-^)}, (1.2) 

T 

fs = - lim —\nZN,s, (1-3) 

m. = - — . (1.4) 

A linear susceptibility Xi,5 and the third order susceptibility X3,5 are defined as follows: 
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X3,S 






(1.5) 
/i=0 



(1.6) 

h=0 



3! a/i4 

In §2 we calculate a scaling function of the magnetization in the classical limit. Moreover 
we extend this scaling function to the case of an arbitrary spin. In §3 to check our extended 
scaling function for S'=l/2, we calculate the scaling function numerically by the Bethe 
Ansatz method. We calculate the Helmholtz free energy and the third order nonlinear 
susceptibility. We compare numerical results with the extended scaling function. In §4 we 
summarize the scaling function of the spin S ferromagnetic Heisenberg chain and new results 
of the numerical calculation for 5*= 1/2. In Appendix we calculate the third order nonlinear 
susceptibility for the classical Heisenberg chain analytically. 

We define the classical Heisenberg Hamiltonian liP with classical spins r^ as follows: 

N N 

H^ = -JoE^^-^m-/^E^- (1-7) 

In the classical chain, a partition function Z^, the Gibbs free energy per site /*", a magnetiza- 
tion m*" and a susceptibilities Xi, Xs are defined by changing from H to H^ in eqs.( |1.2| )-( |1.6| ). 



II. THE SCALING FUNCTION AT AN ARBITRARY SPIN 

A. Analytical Calculations in the Classical Limit {S — > cxo) 

Let's consider the classical limit (i.e.S* -—>■ oo) of the ferromagnetic Heisenberg chain in a 
magnetic field. We define normalized spin operators Sj as follows: 



5, 
S' 



(2.1) 



By using Sj's, the Hamiltonian defined by eg . (11.11) is written as follows: 

TV N 

H = -JoE^^*m-/iE<' (2-2) 

i=l 1=1 

where 

Jo = J^^ (2.3) 

We can regard the spin operators Sj as classical three dimensional unit vectors r^ in the 
classical limit, because the spin operators Sj become commutable. ( For example, [s^, s^] = 
ihs^/S ^ in the limit of S* ^ oo on condition that JS^ is fixed. ) Therefore, to calculate 
susceptibilities in the classical limit [jlO|, we can use the classical Hamiltonian Ti'" defined by 



eq. (|1.7|) instead of the quantum Hamiltonian of eq.( pl2D . As T/Jq -^ 0, we get susceptibilities 
as follows: 

2 7 

Xi,5=oo -^ :^, (2.4) 

44J3 

They are calculated from the 2-point and 4-point correlation functions at h=0. Details 
are given in Appendix. From eqs. (|2.4|) and (|2.5|) , we get a scaling relation between the 
magnetization, the magnetic field and the temperature as follows: 

ms=oo{T,h)=Fl^:^y (2.6) 

where 

2 44 

F(x) = -x x^ + 0(x^). (2.7) 

^ ^ 3 135 ^ ^ ^ ' 



B. The Scaling Function at an Arbitrary Spin from the Classical Limit 

In the previous subsection we get the scahng function in the infinite spin hmit where 
we can use the classical Hamiltonian to calculate the scaling function. In this subsection 
we consider the finite spin case. At low temperatures, the correlation length becomes large 
and the system is nearly ordered. Excitations behaves like spin waves independently of S. 
This is the reason why the finite S Heisenberg chain belongs to the same universality class 
as the infinite spin chain (i.e. the classical chain). Therefore the classical Heisenberg chain 
becomes a good model to calculate the scaling function of the magnetization for the spin 
S quantum chain. On the basis of this physical picture, we expect the scaling relation of 
the arbitrary spin S Heisenberg chain at sufficient low temperatures and {S'^Jh/T'^ ^ 1) as 
follows: 

ms{T,h) = F(^^y (2.8) 

where the scaling function F{x) is eq.( p.7|) for the arbitrary spin 5*. In the next section we 
check the extended scaling function for S =1/2 by comparing with the numerical result 
given by the Bethe Ansatz method. 

III. THE SCALING FUNCTION AT S = 1/2 

A. Methods of Numerical Calculations at 5 = 1/2 

We calculate numerically the Gibbs free energy and the magnetization of the 5=1/2 
Heisenberg chain in a magnetic field. One of the authors [§] has given the following equations 
with infinite numbers p;'s 



AT 



hi /, In 1 ^^ 

h 7r(/ ) H Xi 

T ^ 2' 2i 

Jpii 



(3.1) 



2Tipf + 1) 



/ = 1, 2, 3, ■ ■ ■ , oo. 



+ J2^nL{pi,pj)L{pi,-pj), (3.2) 



where 



L{x,y) = 



W + 
-iy + 



1 

_ 1 

{l+ix)_ 



(3.3) 



These p^'s give the maximum eigenvalue A°(T/ J, h) of the transfer matrix and the Gibbs 
free energy per site as follows: 



A%T/J,h) = 2l[[- 



■^ ^'-'2 , ^\^^2 , im1/2 



tt] M + mi+w. 



fs=i/2 = f{T/J,h) = -TlnA\T/J,h). 



(3.4) 
(3.5) 



The magnetization is given by a differentiation of the free energy with respect to h. We 
need to calculate dpi/dh. We define g;'s as follows, 



_dpi _ _ dpi 



Differentiating eqs.( p.l[) -( ^73|) , we obtain the following equation for g/'s: 



4T 



T AT{pf + lf 2 [ 



where 



Yi = Yl [(^{PhVMj - P{PhPj)(li - (^(Pi^ -pMj - (^{Ph -Pj)<li] 



a{x,y) 



P{x,y) 



1 



W + {1- T^) -W 



1 



1 ^ ' 



l+ix ' 



(1-ixY 



iy + {l 



1 



+ 



(l+ixY 



) -w + (1 



1 



l—ix^ ^ ^ 1+ix ' 

The magnetization per site is represented by pi's and g^'s: 



7715=1/2 = m = 2Tj2 ^^— 



1=1 Pf + 1 
In actual numerical calculations we can approximate p/'s and g^'s as follows 

Pi = -j-{l ^2^+ «Im(pL), 
qi = Ql, 



(3.6) 
(3.7) 

(3.8) 

(3.9) 
(3.10) 

(3.11) 

(3.12) 
(3.13) 



for / > L, where L is an integer. We can obtain solutions of Pi, ■ ■ ■ ,Pl and gi, ■ ■ ■ , g^. By 
increasing L, we can estimate the product in eq. (|3.4|) and the summation in eq.( p.ll|) . 
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B. Numerical Results at 5 = 1/2 

1. The Expansion of the Free Energy by the Least-Squares Method 

Let us consider the Helmholtz free energy g{t,m) per site of the Heisenberg chain in a 
magnetic field which is defined as follows: 

g{t,m) = f{t,h) + mh, (3.14) 

where 

t^j. (3.15) 

In the case oim = 0, at low temperatures it has been already obtained by several authors 
HJTTH . Takahashi and Yamada [§] obtained 



g{t, 0) = Jt^-^ {-1.042 + l.OOt"-^ - 0.9t + 0{t^-^)] . (3.16) 

Xi,s=i/2 = j^ {0.1667 + 0.581t°-^ + 0.68t + 0(t^-^)} . (3.17) 

From the modified spin wave theory |Tl| the susceptibihty equals: 



Xi,s=i/2 = j^ {(1/6) + 0.582597t°-^ + 0.678839t + 0{t^-^)] . (3.18) 

When m is finite, the coefficients of the free energy and the magnetization will be changed 
from the values of the m=0 case. By analyzing our numerical data, we can get the m 
dependence of the free energy as follows: 
g{t,m) _ ^ ^^1.5 2 



J 



ao{t)t^-^ + r{ai(t)m^ + a2(t)m^ + 0{m^)}, (3.19) 



where 



ao(t) = -1.04218 + l.OOt"-^ - 0.94t + 0.9t' ^ (3.20) 

^U^) ^ o 1 7T7«VT7rc^o«7n"5~r7r«^?o^~TTT77r5Y r ' (3.21) 



2 [ (1/6) + 0.5826t0-5 + 0.678t + 0{t^ 



"^^ '* 4I 0.153 + 0.83t0-5 + 2.4t + 0(ti-5)j' ^' ^ 



We show the t°-^ dependence of g{t,m)/{Jt^^'^) at m = 0, 0.28, 0.632 in Fig.|l|. This figure 
gives us an important result that the m dependence of g{t^m) does not begin with t^'^ but 
t^. We plot ao(t), ai{t) and 02 (t) versus t°'^ in Figs.^^, where solid curves are given by the 
least-squares method and where filled circles denote numerical data. 

Since h = dg/dm, eq. (|3.19|) gives us the relation between m and h at fixed values of t as 



follows: 

h(i mi 

2ai(t)m + 4a2(t)m^ + 0(m^). (3.23) 



When h/{Jt^) ^ 1, we can transform eq.( p.23| ) into the following equation: 



By differentiating eq.( p.24|) by /i, we have the linear susceptibility as follows: 

1 



Xi, 5=1/2 



2Jt'^ai{t) 
= -^ {(1/6) + 0.5826t°-^ + 0.678t + 0(t^-^)} . (3.25) 

This result is consistent with eqs.( p.l7D and (|3.18|) . Moreover we get the third order suscep- 



tibility as follows: 



_ 02 (^) 



JH^ 



— {0.00504 + 0.0431t°-^ + 0.13t + 0{t^- 



(3.26) 



2. The Magnetization Curve 

In eq.( p.24| ), we can get a relation between m and h, at the given temperature t. Here 
we consider how this relation behaves in the limit of t ^ 0. Equation ( |3.11| ) gives us m at 



fixed values of bath the temperature and the magnetic field. In Fig.^, we plot h/{Jt^) versus 
^0.5 ^Qj^ g vaiygg oi m : m = 0.01,0.04,- ■ -, 0.63. From this figure, it is clear that ln(/i/ Jt^) 



scale with t°'^. The extrapolated values given by Fig.^ are shown by empty squares in Fig.|^, 
where a bold curve is drawn as a guide for eyes. Moreover in this figure we show m versus 
ln(/i/Jt^) for 5 values of t : t = 0.2,- ■ ■ ,0.005 by filled circles. This figure tells us that as 
t — > 0, magnetization curves are going to a bold curve, which is regarded as the scaling 
function. Especially for small /i/(Jt^), the scaling relation between the magnetization, the 
magnetic field and the temperature is the following equation: 

. ,x 1 fJh\ (jhV ^fJhV 

^5=i/2(r,/.) = - — -0.00504 — +0 — , (3.27) 



from eqs. (|3.25|) and (|3.26|) . 



C. Comparison ^vith the Extended Scaling Function 



In the previous subsection, we get the scaling function ( |3.27| ) of the S = 1/2 quantum 
Heisenberg chain by the Bethe Ansatz method. In this subsection we compare it with the 
extended scaling function of eq.(|2T 



In the case of S'=l/2 in eq.(p.8|), the extended scaling function is as follows: 
^( Jh\ l(Jh\ 11 (JhV ^fJhV 

The third order coefficient of the scaling function given by the numerical calculation is 
-0.00504 from eq.(^^. Since (11/2160) ~ 0.00509, the extended scaling function of 



eq.( |3.28|) agrees with the numerical results of eq.( |3.27|) within the limits of an error 



IV. SUMMARY AND DISCUSSION 

We propose the scaling relation between the magnetization, the magnetic field and the 
temperature for the spin 5* ferromagnetic Heisenberg chain in the analogy with the classical 
limit as follows: 

ms{T,h) = F(^], (4.1) 



where 

2 44 

F{x) = -x-—x' + 0{x'). (4.2) 

In the case of S=l/2, we calculate the Helmholtz free energy of the S=l/2 ferromag- 
netic Heisenberg chain in a magnetic field by the Bethe Ansatz method. According to this 
calculation the Helmholtz free energy is expanded in respect of both t°'^ and m'^ at the 
low temperatures and the low magnetization. The coefficient of the first term t^'^, which 
is independent of m, agrees with the result given by integral equations P]. We get the m 
dependence of the coefficients of t^^t^'^ and t^ up to the forth power m. This results give 
us the third order nonlinear susceptibility. Moreover we get not only m-h curve but also 
the shape of the scaling function. Especially for small h/{Jt^), we get the coefficients of the 
scaling function up to the third order. We show that this scaling function agrees with the 
proposed scaling function. According to our physical picture we will be able to obtain the 
higer order scaling function analytically of the spin 5* ferromagnetic Heisenberg chain by 
calculating the higer order susceptibilities in the classical limit. 
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APPENDIX A: THE NONLINEAR SUSCEPTIBILITY OF THE CLASSICAL 

HEISENBERG CHAIN 

We consider the classical open Heisenberg chain with classical spins rj which are three 
dimensional unit vectors. We choose a magnetic field h {\h\ = h) as 2;- axis in this three 
dimensional space. The Hamiltonian of the classical Heisenberg chain is defined as follows: 

n^ = -jon^ -hn^, (i.i) 



where 



N 

K = E^^-^m, (1-2) 

i=l 

N 

^f = E^f- (1-3) 



The partition function Z^ is 



z-^l---l^---^eMKn^ + Lnf), (1.4) 



where 



K^^,L^^, (1.5) 

and dQi is the element of solid angle for the vector rj. Then the Gibbs free energy per site 

f = -^\nZ^. (1.6) 



The linear susceptibility x? ^-i^d the nonlinear susceptibility Xs are defined as follows: 

(1-7) 



X? 



Q2fC 



xi 



dh^ 



1 d^f^ 



3! dh^ 



h=0 

(11 



h=0 

To evaluate these susceptibilities, we must calculate the 2-point and 4-point functions of 
this classical Heisenberg chain without a magnetic field {i.e.h = 0), which are defined by 

11 



* J 



< r^r^rlrf > 



Tirfr] exp{Kn^) 
TiexpiKH^) ' 

TiexpiKH^) 



By using eqs.(p..4|)-( |1.10| ), the susceptibilities are derived as follows: 



X? 



NT 



A{K), 



(1.9) 

:i.io) 



(1.11) 



X3 



3\NT^ 



BiK), 



where 



a{i,j,k,l) 



N N 

N N N N 

i=ij=i fc=i 1=1 



< r-r^rlrf > - < r-r^ >< r^rf > 



- < r-rl >< r^rf > - < r-rf >< r^-rl > 



Fisher [ffUi gave the 2-point function as follows: 



< r^r^j >-- 



:i.l2) 



:i.i3) 
:i.i4) 



:i.i5) 



;i.i6) 



where u is Langevin function and u = cothi^ — {1/K). Then we have A{K) as follows: 

A{K): 



u 



J-« 



2E + E 
As T/Jq -^ (i.e. K -^ oo), Langevin function behaves as follows: 



u^l--. 



;i.i7) 



:i.i8) 



Therefore eqs.( |l.ll| ),( |1.17| ) and ( p.. IS] ) give us the hnear susceptibility []TU[ at low tempera- 
tures as follows: 
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X? 



2J^ 
3T2' 



:i.i9) 



Next we calculate the 4-point correlation function for the classical Heisenberg chain 
without a magnetic field. Let the polar angle of a spin rj referred to the z-axis be Bj. We 
define 6'j_|_i and 0j+i as polar coordinates for r^+i ,which are referred to Vi as the polar axis 
with ri defining the reference plane for 0i+i. We have 



cos Gj+i = cos 9j cos 6'i+i + sin 6j sin 6*4+1 cos (pi+i- 



;i.2o) 



Here we show two recurrence relations to evaluate the 4-point correlation function. Let i, 
j, and k be such integers as 1 < i < j < A; < A^. By using eq.( |L20| ) and < cos0j >= 0, the 
first recurrence relation is 



< cos 9j • • • cos 9j cos 6fc > = < cos 6j ■ ■ ■ cos Qj cos Qk-i >< cos 9k > 

= < cos 9i ■ ■ ■ cos^ Qj > u^~^ . 



;l2i) 



We square eq.(|1.20|) and we obtain the second recurrence relation as follows: 



< cos^ 9j cos^ 0j > = < cos^ Oj cos^ ©j-i > v + 



u 
3K 



1 (Av^- 



9 



+ 1 



(1.22) 



where v 
follows: 



1 - {3u/K). We get the 4-point function by eqs.dmil), ([T^ ) and ([T^ ) as 



< i^i^jr^rf > = < cos 9j cos Qj cos 6^ cos 6/ > 



u^ ' < cos^ Qi cos^ Qk> u^ ^ 



-u 



j-t 



'4yk-j 



9 \ 5 
By using eqs.( p..l5| ), (|L16|) and (|1.23| ) we have 



+ 1U' 



k 

l-k 



i{i,j, k, I) = -u^-' i-v^-^ - 2M2(^-i) \ u^-k 

u L 



where i < j <k <l. 



;L23) 



1.24) 
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We consider B{K) to calculate Xs- By using eqs. (|1.14]) and (|1.24|) , we sum up a{i,j, k, I) 
for all cases of {i,j, k, 1} as follows: 



B{K) 

24 E +12 E +6 E +8 E + E 

i<j<k<l i<j=k<l i=j<k=l i=j=k<l i=j=k=l 

V /'l+u\^ 2{5u^ + Au + 1) fl + u 



151 — fVl — M 



151-nP \l-u 



a{ij,k,l) 
N + 0{1). 



;i.25) 



As T/Jq — >• (i.e. K — >• oo), v behaves as follows 

3 



v^l 



Then the asymptotic behavior of B{K) is 



K' 



B(K) ^ -—K^N. 
^ ' 45 



^1.26) 



;i.27) 



Using eqs.(|1.5|), ( |1.12| ) and ( |1.27| ), we get the nonlinear susceptibility for large K (i.e. low 
T) as follows: 



Xs 



c , 44J3 



135T6 ■ 



;i.28) 
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Figure Captions 

1. The t°-^ dependence of g{t,m)/{Jt^-^) for m = 0, 0.28, 0.632. As t -> 0, all curves 
approach to the same value —1.04218 (filled square). The curves are a guide for eyes. 



2. When m = in eq.(§]TP|), ao(t) = g{t, 0)/( Jt^-^). We plot ao(t) versus t°-^ A solid curve 
is given by the least-squares method, which is —1.04218 + 1.00t°-^ — 0.94t + 0.9t^-^. 

3. When we expand h/{Jmt^) with rn^ for fixed t, we can get expansion coefficients. We 
plot the first order coefficient 2ai(t) versus t^'^ by the same method as in Fig.|[ A solid 
curve is given by the least-squares method, which is l/{(l/6) + 0.5826t°'^ + 0.678t + 
0.04ti-5}. 

4. We plot the second order coefficient 4a2(t) versus t°'^ by the same method as Fig.|^. A 
solid curve is given by the least-squares method, which is 1/(0.153 + 0.83t°'^ + 2.4t}. 

5. We plot \n{h/{Jt^)} versus t"-^ for m = 0.01, 0.04, 0.06, 0.1, 0.14, 0.28, 0.5, 0.63. A 
quadratic function is used as our fitting function. Results are expressed by solid lines. 
Their y-intercepts (empty circles) are limitation values of \n.{h/{Jt^)} as i ^ 0. 

6. We plot m versus \\i{h/{Jt^)} by filled circles for t = 0.2, 0.1, 0.05, 0.025, 0.005. We 
draw fine curves as a guide for eyes. When t — ^ 0, we plot the limitation values (empty 
squares) given in Fig.^ A bold curve is drawn as a guide for eyes. As we increase a 
number of the limitation values, the bold curve is getting to the scaling function. 
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